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Abstract 



The purpose of this work is to explore the role that random arbitrage opportunities 
play in pricing financial derivatives. We use a non-equilibrium model to set up a stochas- 
tic portfolio, and for the random arbitrage return, we choose a stationary ergodic random 
process rapidly varying in time. We exploit the fact that option price and random arbi- 
trage returns change on different time scales which allows us to develop an asymptotic 
pricing theory involving the central limit theorem for random processes. We restrict our- 
selves to finding pricing bands for options rather than exact prices. The resulting pricing 
bands are shown to be independent of the detailed statistical characteristics of the arbi- 
trage return. We find that the volatility "smile" can also be explained in terms of random 
arbitrage opportunities. 
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1 Introduction 



It is well-known that the classical Black-Scholes formula is consistent with quoted options 
prices if different volatilities are used for different option strikes and maturities [1 . To 
explain this phenomenon, referred to as the volatility " smile" , a variety of models has been 
proposed in the financial literature. These includes, amongst others, stochastic volatility 
models El El, and the Merton jump-diffusion model [5]. Each of these is based on 
the assumption of the absence of arbitrage. However, it is well known that arbitrage 
opportunities always exist in the real world (see [3IZj). Of course, arbitragers ensure that 
the prices of securities do not get out of line with their equilibrium values, and therefore 
virtual arbitrage is always short-lived. One of the purposes of this paper is to explain the 
volatility "smile" phenomenon in terms of random arbitrage opportunities. 

The first attempt to take into account virtual arbitrage in option pricing was made 
by physicists in [HJ CD]. The authors assume that arbitrage returns exist, appearing 
and disappearing over a short time scale. In particular, the return from the Black-Scholes 
portfolio, II = V — SdV/dS, where V is the option price written on an asset S, is not 
equal to the constant risk-free interest rate r. In |§JEI] the authors suggest the equation 
dU/dt = (r + x (t))n, where x (t) is the random arbitrage return that follows an Ornstein- 
Uhlenbeck process. In [TT| IT2*] this idea is reformulated in terms of option pricing with 
stochastic interest rate. The main problem with this approach is that the random interest 
rate is not a tradable security, and therefore the classical hedging can not be applied. This 
difficulty leads to the appearance of an unknown parameter, the market price of risk, in 
the equation for derivative price JUE]- Since this parameter is not available directly 
from financial data, one has to make further assumptions on it. An alternative approach 
for option pricing in an incomplete market is based on risk minimization procedures (see, 
for example, [TJ d UJ EH Ej ) . 

In this paper we follow an approach suggested in [3] where option pricing with stochas- 
tic volatility is considered. Instead of finding the exact equation for option price, we focus 
on the pricing bands for options that account for random arbitrage opportunities. We 
exploit the fact that option price and random arbitrage return change on different time 
scales allowing us to develop an asymptotic pricing theory by using the central limit the- 
orem for random processes JH] • The approach yields pricing bands that are independent 
of the detailed statistical characteristics of the random arbitrage return. 
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2 Model with random arbitrage return 



Consider a model of (S, B, V) market that consists of the stock, S, the bond, B and the 
European option on the stock, V. To take into account random arbitrage opportunities, 
we assume that this market is affected by two sources of uncertainty. The first source 
is the random fluctuations of the return from the stock, described by the conventional 
stochastic differential equation, 

dS 

— = puit + adW, (1) 

where W is a standard Wiener process. The second source of uncertainty is a random 
arbitrage return from the bond described by 

dB 

— = rdt + £{t)dt, (2) 

where r is the risk-free interest rate. Given that there are random arbitrage opportunities, 
we introduce here the random process, £(£), that describes the fluctuations of return 
around rdt. The same mapping to a model with stochastic interest rate is used in P| I1U[ ITT]. 
The characteristic of the present model is that we do not assume that £ (t) obeys a given 
stochastic differential equation. For example, in [HI CHI IHj follows the Ornstein- 
Uhlenbeck process [19 3 . 

It is reasonable to assume that random variations of arbitrage return £(£), are on the 
scale of hours. Let us denote this characteristic time by r arb . This time can be regarded 
as an intermediate one between the time scale of random stock return (infinitely fast 
Brownian motion fluctuations), and the lifetime of the derivative T (several months): 
<< T ar \j «T . In what follows, we exploit this difference in time scales and develop an 
asymptotic pricing theory involving the central limit theorem for random processes. 

Now we are in a position to derive the equation for V. Let us consider the investor 
establishing a zero initial investment position by creating a portfolio II consisting of a long 
position of one bond, B, ^ shares of the stock, S, and a short position in one European 
option, V, with an exercise price, K, and a maturity, T. The value of this portfolio is 

dV 

V = B + -S-V. (3) 

The Black-Scholes dynamic of this portfolio is given by two equations, dU/dt = 0, 
and II = (see PQ). The application of Ito's formula to (J3J) together with ([TJ and J2J) 
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with £(£) = 0, leads to the classical Black-Scholes equation: 



dV a 2 S 2 d 2 V BV T/ n 

1 \-rS rV = 0. (4) 

m 2 as 2 ds K ' 



The natural generalization of dU/dt = is a simple non-equilibrium equation: 

dU IT 



(5) 



9* Tarb ' 

where r ar fe is the characteristic time during which the arbitrage opportunity ceases to exist 
(see 0). 

By using a self-financing condition, dU = dV — ^dS + dB, and Ito's lemma, one can 
derive from (0) and (J2J the equation involving the option V(t, S) and the portfolio value 

n(t,s), 

Note that this equation reduces to (HJ) when II = and ^(t) = 0. 

To deal with the forward problem we introduce the non-dimensional time 

t=^^, 0<r<l, (7) 



and a small parameter 

^~arb 



« 1. 



T 

This parameter plays a very important role in what follows. In the limit e — > 0, the 
stochastic arbitrage return £ becomes a function that is rapidly varying in time. It is 
convenient to use the following notation, £ (|) (see [3]). It follows from (jSj) that the value 
of the portfolio II, decreases to zero like exp(— t/eT). Thus, one can assume that IT = 
in the limit e — > 0. We find from (jUJ) that the associated option price, V s (r, S) , obeys the 
following stochastic PDE 



dV e a 2 S 2 d 2 V £ n 3V e p ,jJc,dV £ ,. , 

+ rS—--rV £ + £(-)(S—--V £ \ (9) 



dr 2 8S 2 dS S V V dS 

subject to the initial condition: 

V e (0,S) = max(S- K,0) (10) 

for a call option, where K is the strike price. Here, for simplicity, we keep the same 
notations for the non-dimensional volatility a and the interest rate r. The same PDE is 
used in but with the Ornstein-Uhlenbeck process for £. 
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3 Asymptotic analysis: pricing bands for the options 



To analyze the stochastic PDE (jHJ), we have to specify the statistical properties of the 
random arbitrage return £. Suppose that £(r) is a stationary ergodic random process 
with zero mean, < £(r) >= 0, such that, 

POO 

D= / <£(r + S )£(r) >ds, (11) 
Jo 

is finite. The key feature of this paper is that we do not assume an explicit equation 
for £(t) unlike the works QUI E], where the random arbitrage return £ follows the 
Ornstein-Uhlenbeck process. 

According to the law of large numbers, V £ (r, S) converges in probability to the Black- 
Scholes price, Vbs {t, S) , as e — > 0. One can split V £ (r,S) into the sum of the Black- 
Scholes price, Vbs ( r ; S), and the random field Z £ (r, S), with the scaling factor y/e, giving 

V* (r, S) = V BS (t, 5) + ^ (r, 5) . (12) 

Substituting (fT2|) into 0, and using the equation for Vbs (t, S), we get the following 
stochastic PDE for Z £ (t, S), 

^-^ + (-«;ii(«l-i + f(^-S' ™ 

Our objective here is to find the asymptotic equation for Z £ (t, S) as e — > 0. One can see 
that Eq. ( llflj) involves two stochastic terms proportional to and £ _1 / 2 £(^). Ergodic 
theory implies that the first term in its integral form converges to zero as e — > 0, while 
the second term converges weakly to a white Gaussian noise rj (r) with the correlation 
function 

<^(ti)77(t 2 ) >= 2D5 (n - r 2 ) , (14) 

where the intensity of white noise, D, is determined by (jlljl (see [TB| and Appendix A). 
Thus, in the limit e — > 0, the random field Z £ (r, 5*), converges weakly to the field Z(t, S) 
that obeys the asymptotic stochastic PDE: 

dZ 1 2c2 d 2 Z , , c dZ „ , (JV BS 



a; = r 25 w + r(5 M - z) + ( s i>f - ^ J " w - (15) 

with the initial condition Z(0, S) = 0. This equation can be solved in terms of the classical 
Black-Scholes Green function, G(S, Si, t, ti), to give 



/t poo / gy \ 

y (?(£, Si, r, n) ( S i-Qg-- ~ Vbs (n, St) J ^(r^dStdn, 



(16) 



G(S, S 1 , r, ri ) = — - 7 ====e (17) 



where 



(see, for example, [22J). 

It follows from (fT7)|) that since 77 (i) is the Gaussian noise with zero mean, Z (r, S) is 
also the Gaussian field with zero mean. The covariance 

R (r, S, Y) =< Z (r, S) Z (r, Y) > (18) 

satisfies the deterministic PDE: 

dR 1 2a2 d 2 R 1 2 ^ 2 <9 2 i? BR nr> , 

57 = r 5 V + ^ + r{s ds + Y w - 2R)+ 
2D { s i^r ~ yBS (r ' 5) ) ( r ~ yBS (r ' Y) ) (19) 

with -R(0, S, Y ) =0 (see Appendix B). The pricing bands for the options for the case of 
arbitrage opportunities can be given by 

V BS (t, S)±2y/eU (r, S), (20) 

where 

U (r, S) = R (r, S, S) . (21) 

The variance U (r, S) quantifies the fluctuations around the classical Black-Scholes price. 
It should be noted that it is independent of the detailed statistical characteristics of the 
arbitrage return. The only parameter we need to estimate is the intensity of noise D which 
is the integral characteristic of random arbitrage return (see (jlljO . One can conclude that 
the investor who employs the arbitrage opportunities band hedging sells the option for 

V BS (T,S) + 2y/eU (t,S), (22) 

where U (r, S) can be found from ()16|) or (fTTJj) to be 

U (r, S) = 2D j\£ G{S, S u r, n ) (si^- - V BS ^j dS 1 fdr l (23) 

(see Appendix C). 

One can see from (JT5J) or (|2*3j) that the large fluctuations of Z(r, S) occur when the 
function (£^j§ £ — Vbs) takes the maximum value. That is, S d2 ^ s = (the Greek 
r = 0). 
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4 Numerical results 



To determine how the random arbitrage opportunities affect the option price, we solve 
equation (JT§j) numerically. Figure 1 shows a graph of the variance, U(r, S) = R (r, S, S), 
as a function of both time, r, and asset price, S. From this graph, we observe that 
the uncertainty regarding the option value is greater for deep-in-the money options. This 
finding is consistent with the empirical work given in (201- in Figure 2, we plot the effective 
option price given by (|2*2*|) for e = 0.1, and compare it with the Black-Scholes price. Note 
that deep-in-the money options deviate the most from the Black-Scholes option price. As 
we move near at-the-money options the deviation decreases. 

Now we are in a position to discuss the "smile" effect, that is, the implied volatility 
is not a constant, but varies with strike price K and time r. Let us denote the implied 
Black-Scholes volatility by a £ (K, r) . The formula f!22D for the effective option price implies 

V BS (r, S; a £ (K, r),K) = V BS (r, S; a, K) + 2y/eU (t,S;K). (24) 

This equation can be solved for a £ (K, r) with S and e fixed. It follows from (|22|) that 
a £ (K, t) — > a as e — > 0. Fig. 3 illustrates the 'smile' effect, when implied volatility o £ (K, r) 
increases for deep-in and out-of-the money options. Similar results are given in [21] where 
the characteristic time T ar & depends on moneyness Numerical results give a similar 
smile curve where the implied volatility becomes greater as we move away from at-the- 
money. 

5 Conclusions 

In this paper, we investigated the implications of random arbitrage return for option 
pricing. We extended previous works by using a stationary ergodic process for modelling 
random arbitrage return. We considered the case where arbitrage return fluctuates on a 
different time scale to that of the option price. This allowed us to use asymptotic analysis 
to find option pricing bands rather than the exact equation for option value. We derived 
the asymptotic equation for the random field that quantifies the fluctuations around the 
classical Black-Scholes price and showed that it is independent of the detailed statistical 
characteristics of the arbitrage return. In particular, we showed that the risk from the 
random arbitrage returns is greater for deep-in and out-of-the money options. This gives 
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an explanation of the 'smile' effect observed in the market in terms of the random arbitrage 
return. 



Appendix A 

The random process x s (t) defined as 

af(r) = \ [ £(-)ds = e5 [' £(s)ds (A-l) 
£2 Jo e Jo 

converges weakly to the Brownian motion B (r) as e — > 0. It means that e -1 ' 2 £(~) con- 
verges weakly to a white Gaussian noise rj (r) with the correlation function < T)(ti)t)(t2) >= 
2D5 (n - r 2 ) (note that V (r) = ^). 

The purpose of this Appendix is to show that 

lim < [x £ (t)} 2 >= 2Dt. (A-2) 

It follows from (lA-ljl that 

a X. ft — 

< [x £ (r)} 2 >=e 1 r < £(si)£(s 2 ) > d Sl ds 2 . (A-3) 



JO Jo 

By using the well-known formula for stationary process £(r) with zero mean 

PT PT PT PT 

/ / < i{si)i{s 2 ) > dsids 2 = 2t I <i(r + s)i{T)> ds-2 I s < f(r + s)£(r) > ds, 
Jo Jo Jo Jo 

(A-4) 

one finds that 

E{[x £ (t)] 2 } = 2r f E < £(r + s)£(r) > ds - 2e I s < f (r + s)£(r) > ds. (A-5) 
In the limit e — > the second term tends to zero, and so we must have 



\imE{x £ (r)} 2 = 2Dr, 



where 



POO 

D= < £(r + S )£( T ) > ds. (A-6) 
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Appendix B 



In this Appendix we derive the equation for the covariance 

R(r,S,Y) =< Z(r,S)Z(r,Y)> . (B-l) 

First, let us find the derivative of R (r, S, Y) with respect to time 

^=<Z (r, S) — ^— + Z (r, Y) — > . (B-2) 

Substitution of the derivatives dZ ^ S ' > and 9Z< £ Y * > from equation (|15|) into ()B-2|) and av- 
eraging give 

di? 1 2a2 d 2 R 1 2 ^ 2 d 2 # 0i2 nn . 

S^f - V BS \ (Z (r, K) ^(t)> + (V^ - Vbs^ (Z (r, 5) ^(t)> . (B-3) 

This equation involves the correlation functions (Z (r, Y) t)(t)) and (Z (t, S) rj(r)) that 
can be found as follows. Since the random process r\ (r) is Gaussian, one can use the 
Furutsu-Novikov formula to find 

(Z (r, Y) rj(r)) = J (v(r)v(n)) x (^fj^p) ^ ( B " 4 ) 

(see [23J). By using delta-correlated function for rj (r) given by (|T1|). and equation ((T3J) , 
we can find the variational derivative 

(see El) and then the correlation function (IB-41) 



(Z (r, F) ^(r)> = -D ( F ^ - V BS (r, Y") ) . (B-6) 

The equation for the covariance R becomes 

dR 1 2a2 d 2 R 1 2 ^ 2 a 2 R , a dR BR nn . 

^ = r s ^ + r Y w- 2+r{s ds +Y w~ 2R) + 

2D (S^f - V BS (r, S)) (Y^ - V BS (r, Y)\ . (B-7) 
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Appendix C 



In this Appendix we briefly discuss the derivation of (J2BJ). Using (jlfijl . one can average 
Z 2 (t, S) as follows 

U(t, S) = (Z 2 (t, S))= / G(S, S u r, n)G(S, S 2 , r, r 2 ) (C-l) 

Jo Jo Jo Jo 

<jS ' L ' S Si - V B s (ti, Si)) (^r-Si - V B s (r 2 , S2) ] (?7 r]( T 2))dS 1 dS 2 dT 1 dT 2 , 



ds x v ' 7 V 95 

where (77 (ti) 77(t 2 )) = 2D5 (ti — r 2 ), and the intensity of the noise .D is determined by 
(lllj) . It follows from the property of the Dirac delta function that 

U (r, S) = 2D J\J™ G(S, S u r, n) f^^i - ^5 (ri, Si)) dS^dn. (C-2) 
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Fig. 1: Variance, U(t,S), with respect to asset price, S, and time r. The option strike 
price K is 20, the volatility o is 0.4, the interest rate r is 0.1, and the constant D is 0.1. 
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Fig. 2: Effective option price (dashed line) and Black-Scholes price (continuous line). Here 
e = 0.1, and constants K, a, r, and D are taken as in Figs. 1. 
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Fig. 3: "Smile" curve: the implied volatility o e (K, r) as a function of the strike price K. 
Here e — 0.1 and constants a, r, and D are taken as in Fig. 1. 
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